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Abstract 

We describe the interplay between d-wave superconductivity and spin density wave (SDW) order 
in a theory of the hole-doped cuprates at hole densities below optimal doping. The theory assumes 
local SDW order, and associated electron and hole pocket Fermi surfaces of charge carriers in 
the normal state. We describe quantum and thermal fluctuations in the orientation of the local 
SDW order, which lead to d-wave superconductivity: we compute the superconducting critical 
temperature and magnetic field in a 'minimal' universal theory. We also describe the back-action 
of the superconductivity on the SDW order, showing that SDW order is more stable in the metal. 
Our results capture key aspects of the phase diagram of Demler et al. (Phys. Rev. Lett. 87, 
067202 (2001)) obtained in a phenomenological quantum theory of competing orders. Finally, we 
propose a finite temperature crossover phase diagram for the cuprates. In the metallic state, these 
are controlled by a 'hidden' quantum critical point near optimal doping involving the onset of SDW 
order in a metal. However, the onset of superconductivity results in a decrease in stability of the 
SDW order, and consequently the actual SDW quantum critical point appears at a significantly 
lower doping. All our analysis is placed in the context of recent experimental results. 
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I. INTRODUCTION 



A number of recent experimental observations have the potential to dramatically advance 
our understanding of the enigmatic underdoped regime of the cuprates. In the present paper, 
we will focus in particular on two classes of experiments (although our results will also have 
implications for a number of other experiments): 

• The observation of quantum oscillations in the underdoped region of YBCO.->^i^>^'^ 
The period of the oscillations implies a carrier density of order the density of dopants. 
LeBoeuf et al.^ have claimed that the oscillations are actually due to electron-like 
carriers of charge — e. We will accept this claim here, and show following earlier 
work-i^, that it helps resolve a number of other theoretical puzzles in the underdoped 
regime. 

• Application of a magnetic field to the superconductor induces a quantum phase tran- 
sition at a non-zero critical field, Hgdw, involving the onset of spin density wave 
(SDW) order. This transition was first observed in La.2-x Sr2.Cu04 with x = 0.144 
by Khaykovich et al.-. Chang et alM^ have provided detailed studies of the spin dy- 
namics in the vicinity of -ffsdw, including observation of a gapped spin collective mode 
for H < Hsdw whose gap vanishes as H y -ffgdw Most recently, such observations 
have been extended to YBa2Cu306.45 by Haug et al.—, who obtained evidence for the 
onset of SDW order at if ~ 15 T. These observations were all on systems which do 
not have SDW order at if = 0; they build on the earlier work of Lake et al— who 
observed enhancement of prexisting SDW order at if = by an applied field in 'La,2-x 
Sr^Cu04 with x = 0.10. 

We begin our discussion of these experiments using the phenomenological quantum theory 
of the competition between superconductivity and SDW order.— li^ii^ii^ The phase diagram 
in the work of Demler et al.— is reproduced in Fig. [H The parameter t appears in a Landau 
theory of SDW order and tunes the propensity to SDW order, with SDW order being favored 
with decreasing t. We highlight a number of notable features of this phase diagram: 

A. The upper-critical field above which superconductivity is lost, Hc2, decreases with 
decreasing t. This is consistent with the picture of competing orders, as decreasing t 
enhances the SDW order, which in turn weakens the superconductivity. 

B. The SDW order is more stable in the non-superconducting 'normal' state than in the 
superconductor. In other words, the line CM, indicating the onset of SDW order 
in the normal state, is to the right of the point A where SDW order appears in 
the superconductor at zero field; i.e. tc(0) > tc- Thus inducing superconductivity 
destabihzes the SDW order, again as expected in a model of competing orders. 
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FIG. 1: FVom Ref. Q: Phase diagram of the competition between superconductivity (SC) and 
spin density wave (SDW) order tuned by an apphed magnetic field H, and a Landau parameter t 
controhing the SDW order (the efi'ective action has a term t(p^, where (p is the SDW order). The 
labels identifying Hc2, H^dw, and tc{0) have been added to the original figure^i^ but the figure is 
otherwise unchanged. The dashed line does not indicate any transition or crossover; it is just the 
continuation of the line CM to identify tc(0). A key feature of this phase diagram is that SDW 
order is more stable in the metal than in the superconductor i.e. tc(0) > tc- 



C. An immediate consequence of the feature B is the existence of the Une AM of quantum 
phase transitions within the superconductor, representing i^sdw, where SDW order 
appears with increasing H. As we have discussed above, this prediction of Demler et 
al.^'^ has been verified in a number of experiments. 

A related prediction by Demler et al.^^ that an applied current should enhance the SDW 
order, also appears to have been observed in a recent muon spin relaxation experiment.— 

A glance at Fig. [1] shows that it is natural to placed the quantum oscillation 
experiments^'^'^i^'^'^ in the non-superconducting phase labeled "SDW". Feature B above 
is crucial in this identification: the normal state reached by suppressing superconductivity 
with a field is a regime where SDW order is more stable. The structure of the Fermi surface 
in this normal state can be deduced in the framework of conventional spin-density-wave the- 



ory, and we recall the early results of Refs. |20||2l| in Fig. [2J Recent studies^^i^^ have extended 



these results to incommensurate ordering wavevectors Q, and find that the electron pockets 
(needed to explain the quantum oscillation experiments) remain robust under deviations 
from the commensurate ordering at (vr, vr). The present paper will consider only the case of 
commensurate ordering with Q = (7r,7r), as this avoids considerable additional complexity. 

The above phenomenological theory appears to provide a satisfactory framework for inter- 
preting the experiments highlighted in this paper. However, such a theory cannot ultimately 
be correct. A sign of this is that within its parameter space is a non-superconducting, non- 
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FIG. 2: (Color online) Fermi surface evolution in the SDW theory^Si^i. Panel (d) is the "large 
Fermi surface" state appropriate for the overdoped superconductor. The SDW order parameter, ip, 
desribes ordering at the wavevector Q = (vr,7r), and mixes fermion states whose wavevectors differ 
by Q. This leads to the SDW metal state with electron (red) and hole (blue) pockets in panel (b), 
which is the state used here to explain the quantum oscillation experiments.— '^'^I'^i^'^ 



SDW normal state at i7 = and T = (not shown in Fig. [T]). Indeed, such a state is 
the point of departure for describing the onset of the superconducting and SDW order in 
Ref. [ij. There is no such physically plausible state, and the parameters were chosen so that 
this state does not appear in Fig. [H Furthermore, we would like to extend the theory to 
spectral properties of the electronic excitations probed in numerous other experiments. This 
requires a more microscopic formulation of the theory of competing orders in terms of the 
underlyirig electrons. We shall provide such a theory here, building upon the proposals of 
j25l . Our theory will not have the problematic if = 0, T = "normal" state of 
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the phenomeno logical theory, and so cannot be mapped precisely onto it. Nevertheless, we 
will see that our theory does reproduce the key aspects of Fig. [H We will also use our theory 
to propose a finite temperature phase diagram for the hole-doped cuprates; in particular, 
we will argue that it helps resolve a central puzzle on the location of the quantum critical 
point important for the finite temperature crossovers into the 'strange metal' phase. These 
results appear in Section [IV] and Fig. [101 

The theory of superconductivity^ mediated by exchange of quanta of the SDW order 
parameter, has been successful above optimal doping. However, it does not appear to 
be compatible with the physics of competing orders in the underdoped regime, at least in 
its simplest version. This theory begins with the "large Fermi surface" state in panel (d) of 
Fig. [2], and examines its instability in a BCS/Eliashberg theory due to attraction mediated 
by exchange of (p quanta. An increase in the fiuctuations of (p is therefore connected to an 
increase in the effective attraction, and consequently a strengthening of the superconducting 
order. This is evident from the increase in the critical temperature for superconductivity 
as the SDW ordering transition is approached from the overdoped side (see e.g. Fig. 4 in 
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Ref. |27| ). Thus rather than a competition, this theory yields an effective attraction between 
the SDW and superconducting order parameters. This was also demonstrated in Ref. [ij by a 
microscopic computation in this framework of the coupling between these order parameters. 
It is possible that these difficulties may be circumvented in more complex strong-coupling 
versions of this theory^^, but a simple physical picture of these is lacking. 

As was already discussed in Ref. [IJ, the missing ingredient in the SDW theory of the or- 
dering of the metal is the knowledge of the proximity to the Mott insulator in the underdoped 
compounds. Numerical studies of models in which the strong local repulsion associated with 
Mott insulator is implemented in a mean-field manner do appear to restore aspects of the 
picture of competing orders.— Here, we shall provide a detailed study of the model of the 
underdoped cuprates proposed in Refs. Qj^jsj 25, and show that it is consistent with the 
features A, B, and C of the theory of competing orders noted above, which are essential in 
the interpretation of the experiments. 

As discussed at some length in Ref. [s], the driving force of the superconductivity in the 
underdoped regime is argued to be the pairing of the electron pockets visible in panel (b) of 
Fig. [2j Experimental evidence for this proposal also appeared in the recent photoemission 
experiments of Yang et alr^. In the interests of simplicity, this paper will focus exclusively on 
the electron pockets, and neglect the effects of the hole pockets in Fig. [21 Further discussion 
on the hole pockets, and the reason for their secondary role in superconductivity may be 
found in Refs. HSES- 

The degrees of freedom of the theory are the bosonic spinous {a =t; i), and spinless 
fermions g±. The spinous determine the local orientation of the SDW order via 



where a are the Pauli matrices. The electrons are assumed to form electron and hole pockets 
as indicated in Fig. [2b, but with their components determined in a 'rotating reference frame' 
set by the local orientation of ip. This idea of Fermi surfaces correlated with the local order 
is supported by the recent STM observations of Wise et alr^. Focussing only on the electron 
pocket components, we can write the physical electron operators as^i^ 











Z]-9+ + zl9- 






zi9+ + z*^9- 




zi9+ - z*^9- 



:i.2) 



where Gi = (0, vr) and G2 = (tt, 0) are the anti-nodal points about which the electron 
pockets are centered. We present an alternative derivation of this fundamental relation from 
spin-density- wave theory in Appendix [Al 

Note that when = (1,0), Eq. (11. ip shows that the SDW order is uniformly polarized 
in the z direction with = (0, 0, 1), and from Eq. (11.21) we have C| 



^_^(^giGi r^giG2-r) and 



ci = 9- 



- e' 



■"). Thus, for this SDW state, the ± labels on the 9± are equivalent 
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to the z spin projection, and the spatial dependence is the consequence of the potential 
created by the SDW order, which has opposite signs for the two spin components (as shown 
in Appendix [X]). The expression in Eq. (11 .21) for general is then obtained by performing 
a spacetime-dependent spin rotation, determined by z^^ on this reference state. 

Another crucial feature of Eqs. (11.11) and (11.21) is that the physical observables (p and 
are invariant under the following U(l) gauge transformation of the dynamical variables z^ 
and g±. 

z^^e'^z^ ■ g+-^e-'^g+ ; g.^e'^g^. (1.3) 



Thus the ± label on the g± can also be interpreted as the charge under this gauge transfor- 
mation. This gauge invariance implies that the low energy effective theory will also include 
an emergent U(l) gauge field A^. 

We will carry out most of the computations in this paper using a "minimal model" for 
Za and g± with the imaginary time (r) Lagrangian^*^ 
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(1.4) 



where the fermion action is 



2m* 



fx 



+ 9'- 



1 



2m* 



9+ 
9-1 



;i.5) 



and the spinon action is 

N 



Y^[\[dr-lAr)z^\^+V^ 
a=l ^ 



fV - iA)z^ 




N 



1.6) 



Here the emergent gauge field is A^ = {A^-, A), and, for future convenience, we have general- 
ized to a theory with spin components (the physical case is A^ = 2). The field g imposes a 
fixed length constraint on the z^, and accounts for the self-interactions between the spinous. 

This effective theory omits numerous other couplings involving higher powers or gradi- 
ents of the fields, which have been discussed in some detail in previous work.-i^i^iiS^ It also 
omits the 1/r Coulomb repulsion between the g± fermions-this will be screened by the Fermi 
surface excitations, and is expected to reduce the critical temperature as in the traditional 
strong-coupling theory of superconductivity. For simplicity, we will neglect such effects here, 
as they are not expected to modify our main conclusions on the theory of competing orders. 
Non-perturbative effects of Berry phases are expected to be important in the superconduct- 
ing phase, and were discussed earlier;- they should not be important for the instabilities 
towards superconductivity discussed here. 
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As has been discussed earlier the theory in Eq. (11. 4p has a superconducting ground state 
with a a simple momentum-independent pairing of the g± fermions (g+g^) 7^ 0. Combining 
this pairing amplitude with Eq. (11.21] . it is then easy to see^i^ that the physical fermions 
have the needed (i-wave pairing signature (see Appendix . 

The primary purpose of this paper is to demonstrate that the simple field theory in 
Eq. (II. 4p satisfies the constraints imposed by the framework of the picture of competing 
orders. In particular, we will show that it displays the features A, B, and C listed above. 
Thus, we believe, it offers an attractive and unified framework for understanding a large 
variety of experiments in the underdoped cuprates. We also note that the competing order 
interpretation of Eq. (II. 4p only relies on the general gauge structure of theory, and not 
specifically on the interpretation of g± as electron pockets in the anti-nodal region; thus it 
could also apply in other physical contexts. 

Initially, it might seem that the simplest route to understanding the phase diagram of 
our theory Eq. (11.41) is to use it to compute the effective coupling constants in the phe- 
nomenological theory of Ref. Q- However, such a literal mapping is not possible, because, 
as we discussed earher, the phenomenological theory does have additional unphysical phases. 
Rather, we will show that our theory does satisfy the key requirements of the experimentally 
relevant phase diagram in Fig. [H 

A notable feature of the theory in Eq. ( 11. 4p is that it is characterized by only 2 dimen- 
sionless couplings. We assume the chemical potential is adjusted to obtain the required 
fermion density, which we determine by the value of the Fermi wavevector kp- The effective 
fermion mass m* and the spin-wave velocity then determine our first dimensionless ratio 

ai = (1.7) 

■m*v 

Although we have inserted an explicit factor of h above, we will set ^ = fc^ = 1 in most 
of our analysis. Note that we can also convert this ratio to that of the Fermi energy, 
Ep = h?kp/{2m*) and the energy scale m*f^: 

^'^ (1.8) 



m*v'^ 2 

From the values quoted in the quantum oscillation experiment^ m* = 1.9me and nkp = 
5.1 nm~^, and the spin-wave velocity in the insulator v ^ 670 meV A, we obtain the estimate 
«! ~ 0.76. We will also use 

m*v^ ^ 112 meV (1.9) 

reference energy scale. 
The second dimensionless coupling controls the strength of the fiuctuations of the SDW 
order, which are controlled by the parameter t in Eq. (II. 6p . Tuning this coupling leads to a 
transition from a phase with (za) 7^ to one where the spin rotation symmetry is preserved. 
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We assume that this transition occurs at the value t = tc(0) in the metalhc phase (the 
significance of the argument of will become clear below): this corresponds to the line CM 
in Fig. [TJ Then we can charaterize the deviation from this quantum phase transition by the 
coupling 



Note that 02 < corresponds to the SDW phase in Fig. [H while 02 > corresponds to the 
"Normal" phase of Fig. [H For 02 > 0, we can also characterize this coupling by the value 
of the spinon energy gap A^, in the N = 00 theory, which is (as will become clear below) 



It is worth noting here that our "minimal model" (Eq. (11.41) ) in two spatial dimensions has 
aspects of the universal physics of the Fermi gas at unitarity in three spatial dimensions. The 
latter model has a 'detuning' parameter which tunes the system away from the Feshbach 
resonance; this is the analog of our parameter 02 • The overall energy scale is set in the 
unitary Fermi gas by the Fermi energy; here, instead, we have 2 energy scales, Ep and 

The outline of the remainder of the paper is as follows. In Section [Tll we will consider the 
pairing problem of the g± fermions, induced by exchange of the gauge boson A^. We will 
do this within a conventional Eliashberg framework. Our main result will be a computation 
of the critical field Hc2, which will be shown to be suppressed as SDW order is enhanced 
with decreasing t. Section [Till will consider the feedback of the superconductivity on the 
SDW ordering, where we will find enhanced stability of the SDW order in the metal over 
the superconductor. Section [IV] will summarize our results, and propose a crossover phase 
diagram at non-zero temperatures. 

II. ELIASHBERG THEORY OF PAIRING 

In our mininal model, the charge and spin excitations interact with each other through 
the Afj^ gauge boson. So the gauge fiuctuation is one of the key ingredients in our analysis. 
We begin by computing the gauge propagator, and then we will determine the critical 
temperature and magnetic field within the Eliashberg theory in the following subsections. 

We use the framework of the large N expansion. In the limit N = 00, the gauge field is 
suppressed, and the constraint field g takes a saddle point value {ig = rn?) that makes the 
spinon action extremum in Eq. ( 11. 6p . At leading order, the spinon propagator has the form 




(1.10) 



A, 



(1.11) 



t 



(2.1) 



8 




where k is spatial momentum, Un is the Matsubara frequency. The saddle point equation 
for m is 



47r2 

The solution of this is 



—m a2 + 



du f d'^k 1 



(2.2) 



m = 2Tln 



(2.3) 



which holds for —00 < ^2 < 00. This result is plotted in Fig. [31 Clearly, m is a monotonically 
increasing function of 02 • Recall that the positive ^2 region has no SDW order, and m is 
large here. As we will see below, the value of m plays a significant role in the photon 
propagators. 

The photon propagator is determined from the effective action obtained by integrating 
out the spinous and non-relativistic fermions. Using gauge invariance, we can write down 
the effective action of the gauge field as follows: 



NT^ f d^k 



(2.4) 



As in analogous computation with relativistic fermions in Ref. |30|, we separate the photon 
polarizations into their bosonic and fermionic components: 



D2 = ND2b + D2f. 



(2.5) 
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We use the Coulomb gauge, k ■ A = in the computation. After imposing the gauge 
condition, the propagator of A^- from the above action is 1/Di, while that of Ai is 

We will approximate Dih and by their zero frequency limits. Computation of the 
spinon polarization in this limit, as in Ref. 



301 yields 



T 

^ dx^m^ + v^k'xil - x) coth ( ^""^ + ""^'^'"^^^ ^ ] (2.7) 



T / / Tfl \ 

Duik) = --5l„(2Binh(-) 



and 

n^,(Jc\ = - _ 

Stt Jq + f2pa;(l - a;) \ 2T 

For the fermionic contributions, we include the contribution of the g± fermions with effective 
mass m* and Fermi wavevector kp. Calculation of the fermion compressibility yields 



D2i>{k) = -— / dx^====== coth ^ — ^ ^ (2.8) 



D,fik,en) = 2 j 



d^q (nf^(gq-k/2) - ^F(gq+k/2)) 

47r^ (icn + k-q/m*) 



« -, (2.9) 

TT 



where ri^;' is the Fermi function. For the transverse propagator, we obtain from the compu- 
tation of the fermion current correlations 



D2f{k,en) + ^Dy{k,en) 



kl 2 fd^qf^2 (q- k)^^ (raF(eq-k/2) -?T'F(£q+k/2)) 



27rm* 771*^7 Att'^ \ k"^ J (ie„ + k-q/m 

kpl^nl 

7rk 



(2.10) 



Putting all this together, we have the final form of the propagators. The propagator of Ar 
is ^ 

NDuik)+m*/'K ^'^'^^^ 

while that of A^ is 

_ I /2 12) 

F J iVD25(A;) + A;^|e„|/(7rfc)- ^ ' ' 
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A. Eliashberg equations 



We now address the pairing instability of the g± fermions. Both the longitudinal and 
transverse photons contribute an attractive interaction between the oppositely charged 
fermions, which prefers a simple s-wave pairing. However, we also know that the trans- 
verse photons destroy the fermionic quasiparticles near the Fermi surface, and so have 
a depairing effect. The competition between these effects c an be addressed in the usual 

we can anticipate that 
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Eliashberg framework.— Based upon arguments made in Refs. 
the depairing and pairing effects of the transverse photons exactly cancel each other in the 
low-frequency limits, because of the s-wave pairing. The higher frequency photons yield a 
net pairing contribution below a critical temperature which we compute below. 

Closely related computations have been carried out by Chubukov and Schmalian^S on a 
generalized model of pairing due to the exchange of a gapless bosonic collective mode; our 
numerical results for below agree well with theirs, where the two computations can be 
mapped onto each other.— 

The Eliashberg approximation starts from writing the fermion Green function using 
Nambu spinor notation. 

t{Un) = iuJn{l - Z{Un))fo + ef^ + (f).{Un)fi (2.13) 

where f are the Pauli matrices in the particle-hole space. Then self-consistency equation is 
constructed by evaluating the self-energy with the above Green function, which yields the 
following equation: 

f (fk' - ~ 

T,{iuJn) = '^z2 -^G{k',uJm)'D{q,k,uj^- ujn) (2.14) 

UJrn 



Note that the first line is a formal expression, with T){q^k,iujm) being a combination of 
the photon propagator and the matrix elements of the vertex. The equations are therefore 
characterized by the coupling A*°*(ci;„); computation of the photon contribution yields the 



explicit expression^"^ 

A*°*(a;„) = \T{u:n) + \L (2.15) 



>'T{uJn) = ^ „ ^ / dk 



kp f^'"^ v^l - {k/2kFY 



27r2m* 7o ND2b{k) + A;j.|cu„|/(7rA;) 

Al = 



m* f^''^ dk 



2n^kF Jo ^1 - {k/2kF] 



2 



_NDih{k) +m*/'!T_ 



(2.16) 



11 



0.45 
0.40 
0.35 
0.30 
0.25 







4 






r 











0.18 
0.16 
0.14 
0.12 
0.10 
0.08 
0.06 
0.04 



0.04 -0.02 0.00 0.02 0.04 



0.04 -0.02 0.00 0.02 0.04 

^2 



FIG. 4: (Color online) The pairing coupling constants associated with the longitudinal (A^) 
and transverse (At(w„)) gauge interactions. The parameter measures the distance from 
the SDW ordering transition in the metal, as defined in Eq. (ll.lOp . The dotted (red), dot- 
dashed (green), dashed (blue), and continuous (black) lines correspond to af/2 = Ep / {m*v'^) = 
0.16, 0.21, 0.26, 0.29. We show At(w„ = SvrT) with T/{m*v^) = 0.016 for the transverse interaction. 
Note that XxiuJn) function is analytic near Q2 ~ in the magnified scale. 



We have divided the total coupling into two pieces based on the different frequency depen- 
dence of the longitudinal and transversal gauge boson propagators. The frequency indepen- 
dent term will need a cutoff for the actual calculation as we will see below. The typical 
behaviors of the dimensionless couplings Xx{uJn), are shown in FiglH 

The longitudinal coupling A/, is around 0.35, and has a significant dependence upon a2, 
which is a measure of the distance from the SDW ordering transition. Note that Xl is larger 
in the SDW-disordered phase (0:2 > 0): this is a consequence of the enhancement of gauge 
fluctuations in this regime. This will be the key to the competing order effect we are looking 
for: gauge fluctuations, and hence superconductivity, is enhanced when the SDW order is 
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suppressed. 

The transverse gauge fluctuations yield the frequency dependent coupling At(co'„). This 
is divergent at low frequencies witb^i^ XxitOn) ~ IdJ^I"^/^. As we noted earlier, this diver- 
gent piece cancels out between the normal and anomalous contributions to the fermion self 
energy.— 1^ We plot the dependence of At(cu„) on the coupling a2 for a fixed Un in Fig. HI 
As was for the case of the longitudinal coupling, the transverse contribution is larger in the 
SDW-disordered phase. 

The full self-consistent Eliashberg equations are obtained by matching the coefficients of 
the Pauli matrices term by term. 



tuJn{l-Z{uJn)) = -7rT\2X'°\uj^-uJn)^=^== (2.17) 

AM = -r5^A-(u;^-^„)^^£^^ (2.18) 

where A(a;„) is the frequency-dependent pairing amplitude. 

Now we can solve the self-consistent equations to determine the boundary of the super- 
conducting phase. Our goal is to look for the critical temperature and magnetic field, and 
we can linearize the equations in A{un) in these cases; in other words we would neglect the 
gap functions in the denominator. 

Z{uJn) = IH V'sgn(e„)AT(t^n - Cn) 



ttT 



|<:n|<|l^n 



Then the solution of the critical temperature of linearized Eliashberg equation is equivalent 
to the condition that the matrix K{un,u;m) first has a positive eigenvalue, where 

K{Un, UJm) = >^T{.^n " ^m) + Al A(u;„ - U^) - ^n,m " ^ " (2.20) 

TXl 

with the soft cutoff function with cutoff Ep 

^(^") = 1 ^ / /F ^2 (2-21) 

1 + Ci[uJnlEFy 

where ci is a constant of order unity. The cutoff Bp is the highest energy scale of the 
electronic structure, so it is not unnatural to set the cutoff with the scale. With this, the 
numerics is well-defined and we plot the resulting critical temperature in Fig. O 

For comparison, we show in Fig. [H] the results for T^. obtained in a model with only the 
transverse interaction associated with XTiij^n)- We can use this T^. to define an effective 
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-0.04 -0.02 0.00 0.02 0.04 

FIG. 5: (Color online) The critical temperature for superconductivity obtained by solution of the 
Eliashberg equations. The lines are for the same parameter values as in Fig. HI The top plot has 
critical temperature scaled with m*v'^, and the bottom is one scaled with Ep. 



transverse coupling, At, by Td Ep = exp(— 1/At). Using Tc/Ep ^ 0.008 for ^2 ~ in 
Fig. [6l we obtain ~ 0.2. This is of the same order as the longitudinal coupling A^, for 
^2 ~ in Fig. m 

Bigger attractive interactions Xplun) and A^ clearly induces a higher critical temperature 
in the SDW-disordered region. Note that this behavior is different from the one of previous 
SDW-mediated superconductivity.— (See the results of Ref . l27| Fig. [4]; near the critical 
region, Tc shows the opposite behavior there.) We have also compared the plots obtained 
by scaling by m*f ^ and Ep. The dependencies on the parameter ai are reversed in two 
plots in the SDW-disordered region. To interpret ai as the doping related parameter, we 
should choose the scaling by m*f ^ because the mass m* and spin wave velocity v are not 
affected much by doping. With this scaling (the first plot in Fig. [5]), the critical temperature 
rises with increase doping at fixed a2] of course, in reality, a2 is also an increasing function 
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FIG. 6: (Color online) As in the top panel of Fig. [5l but with only the transverse pairing interaction, 
AT(wn), included. 
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of doping. 



B. Critical field 

This subsection will extend the above analysis to compute the upper-critical magnetic 
field, Hc2 at T = 0. We will neglect the weak Zeeman coupling of the applied field, and 
assume that it couples only to the orbital motion of the g± fermions. This means that Lg 
in Eq. ( 11. Sp is modified to 



^3 = 9\ 



+ 9l 



{dr - ia;] 

(dr + iAr) 



1 



2m 
1 



2m 



-{y -iA-i{e/c)3.f - li 
V + iA - i{e/c)aif - fi 



9+ 
9-, 



(2.22) 



where V x a = if is the applied magnetic field. 

Generally, the magnetic field induces non-local properties in the Green's function. How- 
ever, in the vanishing gap limit, Helfand and Werthamer proved the non-locality only ap- 
pears as a phase factor (see Ref. l4ll ). The formalism has been developed by Shossmann and 
Schachinger— , and we will follow their method. As they showed, in the resulting equation 
for Hc2, the magnetic field only appears in the modification of the frequency renormalization 

Z{Un). 

The Eliashberg equations in zero magnetic contain a term \unZ{un)\, which comes from 
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the inverse of the Cooperon propagator type at momentum q = 0, C(c<j„, 0), where 
C{uJr,,q) 



2n {-iujnZ{uJn) + E + vfQ cos 6){iuJnZ{uJn) + e) 

27cN{0) 



(2.23) 



where A^(0) is the density of states at the Fermi level per spin. 

Now we discuss the extension of this to if = 0, as described in Refs. I40|j4l|]42l . For this, 
we need to replace C(u;„, 0) by the smallest eigenvalue of the operator 

CM = [d'p [ 0CK,g)e*'i ''e-^''-^ (2.24) 



where tt = p — (2e/7ic) A(f). Using Eq. (22) from Ref. |40|, we find the smallest eigenvalue 
of C{un) is 

CoM = / pdp gdgJo(gp)C(^n,g)e-^'/(2.|,) 
Jo Jo 

POO 

= rl / gc/ge-^'^H/2c(^„,g) (2.25) 







where th = \Jhc/2eH is the magnetic length. 

So the only change in the presence of a field is that the wavefunction renormalization 
Z{un) is replaced by Zu^Un), where 

2|^„|4 / qdq—===^== (2.26) 



ZhM ' Jo \/4u;lZ^M + vlq^ 

= 2\ujn\ / xdx = (2.27) 

We can now insert the modified Z{un) into Eq. (12.271) into Eq. f l2.20p . and so compute 
Hc2 as a function of both ai and the SDW tuning parameter a2- The natural scale for the 
magnetic field is 

= (^^^ kl ^ 534 Tesla, (2.28) 

where in the last step we have used values from the quantum oscillation experiment^ quoted 
in Section [H Our results for Hc2/Hm are shown in Fig. [71 We can see that the critical 
field dependence on a2 is similar to the critical temperature dependence: it is clear that 
SDW competes with superconductivity, and that Hc2 decreases as the SDW ordering is 
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FIG. 7: (Color online) The upper critical field Hc2 as a function of ai and 02 using the same 
conventions as in Fig. HI The magnetic field is measured with the units induced by the fermion 
mass via Hm defined in Eq. p.28p . 



enhanced by decreasing 02- Also, we can compare this with the phenomenological phase 
diagram of FiglH the critical field line in Fig. [7] determines the line B-M-D within Eliashberg 
approximation. Finally, the values of iJc2 in Fig- El are quite compatible with the quantum 
oscillation experiments.-'^i'^i^'^'^ 



III. SHIFT OF SDW ORDERING BY SUPERCONDUCTIVITY 

We are interested in the feedback on the strength of magnetic order due to the onset 
of superconductivity. Rather than using a self-consistent approach, we will address the 
question here systematically in a expansion. 

We will replace the fermion action in Eq. (11. 5p by a theory which has N/2 copies of the 
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electron pockets 



N/2 



a=l 



+ 9-a 



2m* 



2m* 



A g+ag-a - A glagia 



(3.1) 



Here we consider the gauge boson fluctuation more rigorously in the sense of accounting 
for full fermion and boson polarization functions. But we will treat the fermion pairing 
amplitude A as externally given: the previous section described how it could be determined 
in the Eliashberg theory with approximated polarization. 

The large expansion proceeds by integrating out the Za and the g±a, and then expanding 
the effective action for g and - formally this has the same structure as the computation 
in Ref. |30|, generalized here to non-relativistic fermions. At = oo, the g±a and Za remain 
decoupled because the gauge propagator is suppressed by a prefactor of 1/A^. So at this 
level, the magnetic critical point is not affected by the presence of the fermions, and appears 
ai t = where 



1 



duj(Pk 



1 



(3.2) 



87r^ uj^ + 

We are interested in determining the 1/A^ correction to the magnetic quantum critical 
point, which we write as 

note that in the notation of Fig. [H tc = ^c(A). The effect of superconductivity on the 
magnetic order will therefore be determined by -F(A) — -F(O), which is the quantity to be 
computed. The shift of the critical point at this order will be determined by the graphs in 
Fig. 3 of Ref. |30|, which are reproduced here in Fig. [HI Evaluating these graphs we find 



F(A) 



(Pqduj f (Ppde 



■^tt'^ (e^ + t>2p2)2 



(2e + u) 



+ 



+ 



Di{q,iu) \{e + Luy + v'^{p + qy uP- ^ v^(f- 



[D2(g, cu) + (c^Vg2)Di(g, u)\ V (e + ujf + t;2(p + g)2 

1/1 1 



ng(g,c<j) \ijj'^^v^(f (c<j + e)^ + f ^(p + g)^ 



(3.4) 



where l/ng(g, tu) = uS^ + v^q^ is the propagator of the Lagrange multiplier field q. The 
last term involving 11^ is independent of A, and so will drop out of our final expressions 
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(a) 




(c) (d) ^ 




(f) 




FIG. 8: Feynman diagrams for the self energy of Za from Ref. 3Qi. The full line represents z^, the 
wavy line is the propagator, and the dashed line is the g propagator which imposes the length 
constraint on z^. 

measuring the influence of superconductivity: we will therefore omit this term in subsequent 
expressions for -F(A). 

It is now possible to evaluate the integrals over p and e analytically. This is done by using 
a relativistic method in 3 spacetime dimensions. Using a 3-momentum notation in which 

= {vpi,e) and = {yqi^uj) and fp = v^f decPp / (Sn^) , some useful integrals obtained 
by dimensional regularization are: 

1 



p4 











3 



pP^{p + Qy leg 



+ (3-5) 



pP\P + Qf 32Q 32Q 

While some of the integrals above appear infrared divergent, there are no infrared divergen- 
cies in the complete original expression in Eq. (13. 4p . and we have verified that dimensional 
regularization does indeed lead to the correct answer obtained from a more explicit subtrac- 
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tion of the infrared singularities. Using these integrals, we obtain from Eq. (13.41) 



F(A) 



(Pqduj 



^7r3 8(cu2 + t;2g2)i/2 



1 



+ 



cu2 + t;V)^i(g,^) g2/^2(g,cu)+cu2Di(g,cu) 



.(3.6) 



The above expression was obtained in the Coulomb gauge, but we have verified that it is 
indeed gauge invariant. 

We can now characterize the shift of the critical point in the superconductor by determin- 
ing the spinon gap, A^, at the coupling t = tc(0) where there is onset of magnetic order in 
the metal i.e. the spinon gap in the superconductor at H = at the value of t corresponding 
to the line CM in Fig. [H To leading order in this is given by 



1 



1 



m*JM 



(3.7) 



This expression encapsulates our main result on the backaction of the superconductivity of 
the g± fermions, with pairing gap A, on the position of the SDW ordering transition. 

Before we can evaluate Eq. (13.71) . we need the gauge field propagators -Di,2- For com- 
pleteness, we give explicit expressions for the boson and fermionic contributions by writing 



D2{q,uj) + ^Di{q,iu) = D^q^u) + D^^{q,uj) 
D,{q,uj) = Di{q,uj) + Di{q,u). 



(3.8) 
(3.9) 



We can read off the bosonic polarization functions D^i^rp[q^ui) from the exact relativistic 
result of Ref. 0, and the Eq. (El). 
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1 



q^ 



D!^{q,cu) 

^^^^'""^ 16 ^^V + ^' 
For the fermion contribution, let us introduce the Nambu spinor Green's function 



(3.10) 
(3.11) 



1 / iuJ + iq -A 



71 n — iuj 



lm[g{q,n)] 



lm[g{q,n)] 



2E. 



-A 



-A Q-^, 



(3.12) 
(3.13) 
(3.14) 
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where = q'^/{2m*) — fi and Eg = ^J^^ + A^. With the matrix elements of longitudinal 
and transverse parts, the polarizations of the fermions are as: 



D{iq,uj) = - J -^^^-^tr [g{k,e)g{q + k,uj + e)] 
(Pk [ dVL' dVLnpiVL') -npiVL) 



(27r)2 J TT TT iuj + n-Q' 



tr \hnxj{k, VL)luig{q + k, Q!)] 



d^k 1 / _ e^a+g + A^ X / 2{Ek + Ek+q) . 



(27r)2 2V ^fc^fe+, ;V^' + (^fc + ^feV 



r^/- /" d'^k 1 A 2 ■ qY\ f dt , 

D^para = (^(^^^^^ " J 7 j^^tT ^k, e)rsg{q + k, u + e)n] 

d^k khin^e f dQdQ' npirt') -npifl) 



(27r)2 (m*)2 J 71 71 iio + n-n' 

X tr [Im(7(A;, fi)f3lm(7(g + A;, i7')'r3] 
d^k khin'e i / _ aa+g-A^ \ / 2{Ek + Ek+q) . ^2 



(27r)2 (m*)2 2 V ^fc^fc+g J W + (Ek + Eu+^f 

= 5, (3.17) 

where p/ is the density of the fermions and is a Pauli matrix in the Nambu particle- hole 
space. With these results we are now ready to evaluate Eq. (13 .7^ . 

One of the key features of the theory of competing orders was the enhanced stability of 
SDW ordering in the metallic phase. This corresponds to feature B discussed in Section [H 
tc{0) > tc in Fig. [TJ In the notation of our key result in Eq. (13. 7p . where tc(A) = tc, this 
requires A^, > 0. We show numerical evaluations of Eq. fl3.7l) in Fig [3 and find this indeed 
the case. (The values of A used in Fig. [9] are similar to those obtained in Section [TTl near the 
SDW ordering critical point.) Indeed, the sign of A^ is easily understood. In the metallic 
phase, the gauge fluctuations are quenched by excitations of the Fermi surface. On the other 
hand, in the superconducting state, this effect is no longer present: gauge fluctuations are 
enhanced and hence SDW ordering is suppressed. Note that the fact that the g± fermions 
have opposite gauge charges is crucial to this conclusion. The ordinary Coulomb interaction, 
under which the g± have the same charge, continues to be screened in the superconductor. 
In contrast, a gauge force which couples with opposite charges has its polarizability strongly 
suppressed in the superconductor, much like the response of a BCS superconductor to a 
Zeeman field. 
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FIG. 9: The energy /S.^ in Eq. (j3.7p determining the value of the shift in the SDW ordering critical 
point, tc(0) — tc{^). The horizontal axis is the externally given superconducting gap. For numerics 
we fix the parameter ai/2 = Ep/m*v^ = 0.3 

IV. CONCLUSIONS 

This paper has described the phase diagram of a simple 'minimal model' of the under- 
doped, hole-doped cuprates contained in Eqs. (11. 4p . (11. 5p . and (II. 6p . This theory describes 
bosonic neutral spinous Za and spinless charge — e fermion g± coupled via a U(l) gauge field 
A^j_. We have shown that the theory reproduces key aspects of a phenomenological phase 
diagran>i^ii^ of the competition between SDW order and superconductivity in Fig. [T] in an 
applied magnetic field, H. This phase diagram has successfully predicted a number of recent 
experiments, as was discussed in Section [H 

In particular, in Section [Tll we showed that the minimal model had a H(.2 which decreased 
as the SDW ordering was enhanced by decreasing the coupling t in Eq. (11.61) . 

Next, in Section [XTTt we showed that the onset of SDW ordering in the normal state with 
H > Hc2 occurred at a value t = tc{0) which was distinct from the value t = tc{A) in 
the superconducting state with H = 0. As expected from the competing order picture in 
Fig. [H we found tc(0) > tc{A). The enhanced stability of SDW ordering in the metal was 
a consequence of the suppression of gauge fluctuations by the g± Fermi surfaces. These 
Fermi surfaces are absent in the superconductor, and as a result the gauge fluctuations are 
stronger in the superconductor. 

We conclude this paper by discussing implications of our results for the phase diagram 
at T > 0, and in particular for the pseudogap regime above Tc. In our application of the 
main result in Section UTTt tc(0) > tc(A) we have assumed that the A = state was reached 
by application of a magnetic field. However, this result also applies if A is suppressed by 
thermal fluctuations above Tc. Unlike H, thermal fluctuations will also directly affect the 
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FIG. 10: Proposed finite temperature crossover phase diagram for the cuprates. The labels at 
T = are as in Fig. [TJ the onset of SDW order in the superconductor is at tc = td^), while tc(0) 
is a 'hidden' critical point which can be observed only at H > Hc2 as in Fig. [TJ The computations 
in Section HIH show that tc{0) > tc{A). The full line is the phase transition at Tc representing loss 
of superconductivity. The dashed lines are crossovers in the fluctuations of the SDW order. The 
dotted lines are guides to the eye and do not represent any crossovers. Thus, in the pseudogap 
regime at T > Tc the SDW fluctuations are in the 'renormalized classical— (RC) regime; upon 
lowering temperature, they crossover to the 'quantum critical' (QC) and 'quantum disordered' 
(QD) regime in the superconductor. 

SDW order, in addition to the indirect effect through suppression of superconductivity. In 
particular in two spatial dimensions there can be no long-range SDW order at any T > 0. 
These considerations lead us to propose the crossover phase diagram in Fig. [10] in the T, t 
plane. We anticipate that tc(0) is near optimal doping. Thus in the underdoped regime above 
Tc, there is local SDW order which is disrupted by classical thermal fluctuations: this is the 
so-called 'renormalized classical— regime of the hidden metallic quantum critical point at 
tc(0). Going below Tc in the underdoped regime, we eventually reach the regime controlled 
by the quantum critical point associated with SDW ordering in the superconductor, which 
is at tc{A). Because tc(A) < tc(0), the SDW order can now be 'quantum disordered' (QD). 
Thus neutron scattering in the superconductor will not display long-range SDW order as 
T — > 0, even though there is a RC regime of SDW order above Tc. This QD region will have 
enhanced charge order correlations^*^'^'^; this charge order can survive as true long-range 
order below Tc, even though the SDW order does not. Thus we see that in our theory the 
underlying competition is between superconductivity and SDW order, while there can be 
substantial charge order in the superconducting phase. 
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Further study of the nature of the quantum critical point at tc{0) in the metal is an 
important direction for further research. In our present formulation in Eq. fll.4p . this point 
is a transition from a conventional metallic SDW state to an 'algebraic charge liquid— in 
the 0(4) universality class.- However, an interesting alternative possibility is a transition 
directly to the large Fermi surface state.— 

Finally, we note that a number of experimental studies^^'^'^'^i^'^i^ have discussed a 
scenario for crossover in the cuprates which is generally consistent with our Fig. [TUl 
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APPENDIX A: FIELD RELATIONS FROM SPIN DENSITY WAVE THEORY 

This appendix will give a derivation of the relation (11. 2p between the physical electron 
operators and the fields g± and Za using spin density wave theory. This will complement 
the derivation obtained from the doped Mott insulator approach in previous works.-*^ 

We begin the quasiparticle Hamitonian which determines the 'large' Fermi surface in the 
overdoped regime 

Ho = -Yl ^iAaCja = J] ^kC^Cka (Al) 
i<j k 

where we choose the dispersion to agree with the measured Fermi surface. In the presence 
of spin density wave order, at wavevector K = (tt, tt), we have an additional term which 
mixes electron states with momentum separated by K 

^^sdw = ■ X] 4,a'?«/3Ck+K,/3 (A2) 

k,a,/3 

where a are the Pauli matrices. 

Now we focus on the electrons which are near the electron pockets. Let us write 

C(0,7r)a = Cia , C(7r,0)a = C2a (A3) 

and £(o,7r) = ^(7r,o) = ^o- Then, for Neel order polarized as = (0, 0, ip) with > 0, the 
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Hamiltonian for these electrons is 



- ^ (4tC2T - cl;C2| + 4|Cit - 4|Ci J (A4) 



We diagonalize this by writing 

Ho + i/sdw = (^0 - {9I9+ + 9-9-) + (^0 + ^) (h\h+ + hlh_) (A5) 

where 



CiT = {9+ + /i+)/v2 

C2T = (^+ -h+)/V2 

cii = (9- + h-)/V2 

C21 = i-g- + h_)/V2 (A6) 

Now the main approximation we make here is to neglect the higher energy h± fermions. 
We obtain the electron operators for a general polarization of the Neel order as in Eq. (11 .ip 
by performing an SU(2) rotation defined by the Za (and dropping the unimportant factor of 
I/V2) 

':y^^{x) ^ (::)^'^^(-) 

where the SU(2) rotation is 

--(::;?)■ 

These results lead immediately to Eq. (11.21) . In the superconducting state, where {g+g-) 7^ 0, 
they yield 

(ciTCii) = {{\z'i\'^ + \zf) g+g-) 

(c2TC2i) = -{{\z^\' + \zf)g+g-), (A9) 
which implies a d-wave pairing signature for the electrons. 
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